Thermopower with broken time-reversal symmetry 
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We show that when non-unitary noise effects are taken into account the thermopower is in gen- 
eral asymmetric under magnetic field reversal, even for non-interacting systems. Our findings are 
illustrated in the example of a three-dot ring structure pierced by an Aharonov-Bohm flux. 

PACS numbers: 05.70.Ln, 72.20.Pa, 05.70.-a 



In power generation and refrigeration by means of ther- 
mal engines, efficiency plays a basic theoretical and prac- 
tical role. The Carnot bound on efficiency lies at the 
foundations of thermodynamics: for a heat engine func- 
tioning between hot and cold reservoirs at temperatures 
Th and Tc, the efficiency 77, defined as the ratio of the 
output power over the heat extracted per unit time from 
the high temperature reservoir, is upper bounded by the 
Carnot efficiency rjc'. r] < rjc — 1 ^ Tc/Th- 

For systems with time-reversal symmetry, thermo- 
electric power generation and refrigeration is governed, 
within linear response, by a single parameter, the dimen- 
sionless figure of merit ZT ~ [aS"^ / k)T, where a is the 
electric conductivity, S is the thermopower (Seebeck co- 
efficient), K, is the thermal conductivity, and T ^ T^ ~ Tc 
is the temperature. The maximum efficiency is given by 
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Thermodynamics only imposes ZT > and the Carnot 
limit is reached when ZT — >■ 00. 

On the other hand, we have recently shown |1[ that 
for systems with broken time-reversal symmetry the ef- 
ficiency depends on two parameters: a "figure of merit" 
and an asymmetry parameter. In contrast to the time- 
symmetric case, the figure of merit is bounded from 
above; yet the Carnot efficiency can be reached at lower 
and lower values of the figure of merit as the asymmetry 
parameter increases. According to the expression for the 
efficiency, large asymmetry of the thermopower can be 
responsible for highly non-trivial effects [l|, and poten- 
tially can be a useful tuning parameter to control thermo- 
electric efficiency of the material. Hence, finding general 
conditions for asymmetry of the thermopower is of gen- 
eral interest both from practical and purely fundamental 
point of view. 

If time-reversal symmetry is broken, e.g. by means of 
a magnetic field B, then one does not expect the See- 
beck coefficient to be in general symmetric with respect 
to the magnetic field. Yet for the particular case of non- 



interacting systems one has S{B) — S{-'B) as a con- 
sequence of the symmetry properties of the scattering 
matrix [3. Even though this constraint does not apply 
when interactions or inelastic scattering are taken into 
account, and even though there are no general results 
imposing the symmetry of the Seebeck coefficient, the 
latter has always been found to be an even function of 
the magnetic field in purely metallic two-terminal meso- 
scopic systems [3|. On the other hand, Andreev interfer- 
ometer experiments [4I and recent theoretical studies in- 
dicate that systems in contact with a superconductor [5| 
or with a heat bath [y] can exhibit non-symmetric ther- 
mopower. However, accurate numerical simulations of 
various models of two-terminal purely Hamiltonian in- 
teracting dynamical systems, which violate time-reversal 
symmetry, such as a two-dimensional anisotropic and in- 
homogeneous system of interacting particles in a perpen- 
dicular magnetic field [7[, systematically failed to find a 
non-symmetric thermopower. S{B) 7^ S(—B). There- 
fore, it is a remains a completely open and interesting 
problem to understand what requirements must be ful- 
filled in order to actually lead to a thermopower which is 
asymmetric in the magnetic field. 

In this Letter we show that the thermopower is in gen- 
eral asymmetric when non-unitary noise is added to the 
system, even though the system is non-interacting. In- 
deed, in the non-interacting case the symmetry of the 
thermopower is a consequence of the unitarity of the 
scattering matrix, which is broken when noise is added. 
A very convenient way to introduce noise is by means 
of a third terminal, whose parameters (temperature and 
chemical potential) are chosen self-consistently so that 
there is no average flux of particles and heat between 
the terminal and the system. In mesoscopic physics, 
such third terminal, or "conceptual probe" is commonly 
used to simulate phase-breaking processes in partially 
coherent quantum transport, since it introduces phase- 
relaxation without energy damping [8|. We also show 
that, as a consequence of the asymmetry of the Seebeck 
coefficient, a weak magnetic field generally improves ei- 
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FIG. 1. Schematic drawing of tiie model. Tiie third (probe) 
reservoir introduces non-unitary noise. 



ther the efficiency of thermoelectric power generation or 
of refrigeration, the efficiencies of the two processes being 
no longer equal when a magnetic field is added. Our find- 
ings are illustrated by the example of a realistic, asym- 
metric three-dot ring structure pierced by an Aharonov- 
Bohm flux. A main advantage of this model is that it 
can be analyzed exactly, without resorting to approxi- 
mations. 

General setup. The model we consider is sketched 
in Fig. [TJ A system is in contact with left (L) and 
right (i?) reservoirs (terminals) at temperatures Tl = 
T + AT, Tji — T (without loss of generality, we as- 
sume Tl > Tr) and chemical potentials fiL = fJ- + Afi, 
jiR ~ fi. Both electric and heat currents flow along the 
horizontal axis. Non-unitary noise effects are simulated 
by means of a third (probe) reservoir (P) at temperature 
Tp = T + ATp and chemical potential fip = /i + A^ip. 
Let Jpk and Jeu denote the particle and energy cur- 
rents from the fc-th reservoir (k = L, R, P) into the sys- 
tem, with the steady-state constraints of charge and en- 
ergy conservation: J^k Jpk = 0, X^fc ^Ek = 0. The sum 
of the entropy production rates at the reservoirs reads 
•5 = Jlki^Ek - fJ'kJpk)/Tk. Within linear response, 
S ^ J ■ X = ^i^i JiXi, where we have defined the 
4-dimensional vectors J and X: 



J = {eJpL, JqL,e.JpP, Jqp), 
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and where the heat currents Jgk = Jsk ^ l^-Jpk and e is 
the electron charge. The equation connecting the fltrxes 
Ji and the thermodynamic forces Xi within linear irre- 
versible thermodynamics is 
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(4) 



where L is a 4 x 4 Onsager matrix. 

The probe reservoir is adjusted in such a way that J3 = 
J4 = 0, that is, the net particle and heat flow from the 
probe into the system vanishes. It is convenient to write 
Eq. (m in the block matrix form 
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where a stands for (1,2) and /3 for (3,4). The self- 
consistency condition Jp ~ [J^^Ji) — implies Xp = 
—Lpp~ LjjaXa, SO that 

Ja = L Xa, L = Laa — Laj3Lj3i3~ Lfja- (6) 

The problem has then been reduced to two coupled 
fluxes: 
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where the reduced 2x2 Onsager matrix matrix L' fulfills 
the Onsager-Casimir relations 



^JB)^LU-B), (z,j = l,2). 
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We would like to draw the reader's attention to the fact 
that the matrix L' is the Onsager matrix for two-terminal 
noisy transport, with noise modeled by means of a self- 
consistent reservoir. In particular, the Seebeck and the 
Peltier coefficients are given by S" = L[2/{eTL[i) and 
n = L'2i/{eLii). The thermopower is asymmetric when 
L'i2(S)^L^i(B), i.e. U^ST. 

A key point is that, since J3 = .Xi = 0, Ji is the charge 
current from left to right reservoir and the heat is ex- 
tracted from (for power generation) or dissipated to (for 
refrigeration) the left (or right) reservoir only. Therefore, 
wc can apply the analysis developed in Rcf. jl|. In partic- 
ular, the efficiency depends on the asymmetry parameter 
X and on the "figure of merit" parameter y: 



L\ 



12 



Tl ■ 

^21 



y = 



T' T' 

dcti' 



(9) 



For power generation (J2 > and output power lo = 
— JiA/i = —JicTXi > 0) the efficiency i] ~ lo/ J2 has a 
maximum value 
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while for refrigeration ( J2 < 0, w < 0) the maximum of 
the efficiency 77^''^ = J2/0J is 






(11) 



Non-interacting systems. Exact calculation of ther- 
mopower and efficiencies is possible for non-interacting 
models by means of the Landauer-Biittiker approach. We 
start from the bilinear Hamiltonian H ~ Hs + Hp -\- He, 
where the different terms correspond, respectively, to 
the nanoscale electronic system, the reservoirs, and the 
reservoir-system coupling. The tight-binding TV-site sys- 
tem Hamiltonian reads 
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where c„ and cj^ are fermionic annihilation and creation 
operators. The reservoirs are modeled as ideal fermi 
gases: Hr = J2k,q^q'4:q'^kq, where cj.^ creates an elec- 
tron in the state q in the fc-th reservoir. The coupling 
(tunneling) Hamiltonian 



He = J2^^'"l'^kqC^k + t*kqCkqc\j 



(13) 
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establishes the contact between site ik and reservoir k [9| . 
The charge and heat currents from the left terminal 
(reservoir) are given by [10| 



Ji 



/OO 
dEY,[TkL 
-CXD 7. 



(£)/i(B)-Tu(£)/i(£)l, (14) 



r.,M.}r-^<>-f,v^f;,r-'"' 



FIG. 2. Schematic drawing of the three-dot model. 



and Lap ^ -^/Sa, since L^a is obtained from Lap after 
substitution of T^p with TpL and in general T^p ^ Tp^. 
The transmission probabilities are given by [2| 



Tp, = Tr[Tp{E)G{E)T,{E)GHE)l 



(20) 



where the broadening matrices Tk are defined in terms 
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dE{E-fXL)Y.l'^kL{E)fL{E)-TLk{E)fk{E)], of the self-energies E^: rfc(^) = i[Efc(£;) - E^(^)] and 
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(15) 
where fk{E) = {exp[(^ - Hk)/kBTk] + l]-^ is the Fermi 
function and T^i is the transmission probability from ter- 
minal / to terminal k. Analogous expressions can be writ- 
ten for J3 and J4, provided the terminal L is substituted 
by P. 

The Onsager coefficients Lij are obtained from the lin- 
ear response expansion of the currents J^. We have 
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where F{E) = ~Tf'{E) = l/4fcB cosh^[(£: - m)/A:bT]. 
Analogous formulas are obtained for L33, L34 = L43, and 
L44, with the P terminal used instead of L. Note that 
for the non-interacting three-teminal model L12 is still an 
even function of the magnetic field, that is, L12 = ^21- 
On the other hand the symmetry of the off-diagonal ma- 
trix elements is broken for the reduced Onsager matrix 
L' . Indeed, reduction (jH]) involves other off-diagonal ma- 
trix elements of i - between 'left' (1,2) and 'probe' (3, 4) 
sectors - which in general are not even functions of an 
applied magnetic field. The block Lap of matrix L is 
given by 
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Tlp{E)F{E), 
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the (retarded) system Green function G{E) = [E — Hs 
Ek^kiE)]-\ 

Aharonov-Bohm interferometer. As an illustrative, re- 
alistic example we consider a three-dot ring structure 
pierced by an Aharonov-Bohm flux, with dot k coupled 
to reservoir k, as sketched in Fig. [5J The system Hamil- 
tonian reads 

Hs = Efc efe4cfe+ 

(iLfl4cLe^'^/3 + tRPC^pCRe'"^/^ + tpLc\cpe'^l^ + H.c), 

(21) 
and the broadening matrices are Ffe = ''jkclck- We ap- 
ply the Landauer-Biittiker approach to this model, nu- 
merically computing the Onsager coefficients following 

As expected, we obtain asymmetric off-diagonal re- 
duced Onsager matrix elements, that is L'12 7^ ^21; ^^ ^^^ 
as the Aharonov-Bohm flux is non- vanishing and there 
is anisotropy in the systems, for instance when e^ 7^ ep. 
Since the thermopower is not symmetric with respect to 
the magnetic field, i.e. Li2(-B) ^ L^^i^B) = L'^iiB), 
then in general the ratio x = ^'12/^21 7^ 1- The asym- 
metry parameter x can be made arbitrarily small when 
L'21 — ?► or arbitrarily large when L'21 — )- 0, see for in- 
stance Fig. [3l 

Remarks. Large asymmetries ipso facto do not im- 
ply large efficiencies. For example, in the case of Fig. [3] 
when X diverges the figure of merit y and the efflciency 
tend to zero. It is however interesting to compare the 
efficiencies of power generation and refrigeration. While 
in the time-symmetric case the two efficiencies coincide, 
?7max(</' = 0) = 77max('/> = 0), this is uo louger the case 
when X ^ 1. For small fields x is in general a linear 
function of the field, x{(t)) = 1 + q0 + 0(0^), while 
y is by construction an even function of the field, so 
that y{cj)) = 2/(0) + Pcj)"^ + 0(0**). From Eqs. ^ and 
((ni) wc obtain riraa,^{(j)) = ?/max(0)[l + a(j) + 0(0^))] and 
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the Seebeck coefficient with respect to the magnetic field 
afiows in principle, in the linear response regime, to ob- 
tain a finite power at Carnot efficiency. Whether this 
is actually the case remains an interesting open prob- 
lem. An additional interesting open problem is whether 
noiseless interacting systems might exhibit asymmetric 
thermopower. 
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FIG. 3. Ratio x of the off-diagonal matrix elements L12 and 
1/21 of the reduced Onsager matrix at T = 1, ^ = 0.3, el = 0, 
eR — 0.5, ep = 1, all hopping terms tpq = —1, broadenings 
7fe =0.1 independently of energy (wide-band limit). Here- 
after we set e = ft = fes = 1. 



'n'£L^{4') = 77max(0)[l -a(j) + O(02))]. Therefore, a smah 
external magnetic field either improves power generation 
and worsens refrigeration or vice-versa, while the aver- 
age efficiency 77 = [77max(0) + »7mlx (</>)] /2 = ?7max(0) up 
to second order corrections. Due to the Onsager-Casimir 
relations a;(— 0) ~ l/x((p) and therefore by inverting the 
direction of the magnetic field one can improve either 
power generation or refrigeration. 

In conclusion, we have shown that non-unitary noise 
generally leads to a thermopower which is a non- 
symmetric function of the magnetic field. Such general 
result has been illustrated by means of a realistic three- 
dot Aharonov-Bohm interferometer model, which ap- 
pears suitable for experimental investigations by means 
of three-terminal mesoscopic devices. The asymmetry of 
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